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Abstract An outline of Morita's equilibrium ensemble approach to disordered systems 
is given, and hitherto unnoticed relations to other, more conventional approaches in the 
theory of disordered systems are pointed out. It is demonstrated to constitute a general- 
ization of the idea of grand ensembles and to be intimately related also to conventional 
low-concentration expansions as well as to perturbation expansions about ordered refer- 
ence systems. Moreover, we draw attention to the variational content of the equilibrium 
ensemble formulation. A number of exact results are presented, among them general so- 
lutions for site- and bond- diluted systems in one dimension, both for uncorrelated, and 
for correlated disorder. 



1 Introduction 

The present contribution is concerned with the statistical mechanics of systems with 
quenched randomnes. According to a seminal paper by Brout Q, the relevant thermo- 
dynamic potential for the description of such systems is the average of the free energy of 
the system over the distribution describing the disorder configurations k, 

f q = -{PNr l {\nZ N {K)) q . (1) 

Here N denotes the system size, and Z^{n) is the partition function of a system of size 
N at fixed disorder configuration n. 

The evaluation of the quenched free energy (||) is difficult for at least two reasons. First, 
in virtually all configurations k contributing to the average, the system is not translation- 
ally invariant — a circumstance invalidating all conventional calculational tools relying 
on such invariance. Moreover, averaging the logarithm of the partition function usually 
precludes any useful factorization of the averaging process even in situations where such 
factorization is possible when considering the average of Zn(k) itself. 
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For the reasons just outlined only few exactly solved models of systems with quenched 
disorder are known. Where solutions of such models have been obtained, success has as 
a rule been due to simplifying features, such as one-dimensionality of the models (see || 
for a recent overwiew) or of the disorder type || , analytic structures available in spherical 
models [[| or infinite dimensionality where mean-field approximations are exact (see, e.g. 

!)• 

Discomforting to a perhaps even greater degree is the fact that conventional asymptotic 
methods for the study of phase transitions and critical phenomena, such as series expan- 
sions, renormalization group or Monte-Carlo calculations, have in several cases procuced 
conflicting and sometimes even spurious results, when applied to disordered systems. This 
being so, alternative approaches are clearly welcome. 

The present contribution is devoted to one such alternative, which is at the same time 
venerable, not well known and still, we believe, promising. The method is originally due 
to Morita H. The main idea is to treat the configurational degrees of freedom k on the 
same footing as the dynamical variables proper, and to supplement the Hamiltonian of the 
disordered system by a disorder potential <t>(n) which is determined such that configuration 
averaging as implied by (|l]) becomes part of the Gibbs average in an enlarged phase space. 
By this device, it is hoped, the full machinery of equilibrium statistical physics is made 
available for the study of systems with quenched randomness. The idea of an equivalent 
equilibrium ensemble naturally lends itself to the formulation of systematic approximative 
schemes. The perhaps most natural was already suggested by Morita ||; it starts out from 
the so-called annealed approximation, which is then improved by imposing an increasing 
number of constraints on the thermal motion of the configurational degrees of freedom. 
Thus, as pointed out by Falk j7||, annealed approximations, such as employed by Thorpe 
and Beeman || to determine phase boundaries in bond-disordered 2-D Ising models, 
constitute just the first step wihtin Morita's restricted annealing scheme. 

The Morita approach is natural enough to have been used, reconsidered, or rediscovered 
several times. Thorpe || applied essentially this idea to improve upon annealed phase 
boundaries obtained earlier P], a topic taken up again somewhat later by George et al. 



[10], who fell on virtually the same set of ideas, arguing from a somewhat different point of 



view, though. Restricted annealing has been discussed in connection with spin-glasses by 



Toulouse and Vannimenus [11], who also drew attention to formal similarities with lattice 
gauge theories. 



Explicit reference to Morita's original paper is made in |l2j — [17]. Sobotta and Wagner 



[12], restricting themselves to bimodal disorder (see below), streamlined the formal pre- 



sentation of the approach, and applied it in conjunction with the renormalization group 



and the e-expansion to study spin-diluted ferromagnets near dimension 4 [13]. Huber [14] 
discusses the thermodynamic content of the Morita scheme and derives a normalization 
condition for the disorder potential 4>(k) from thermodynamic consistency-considerations. 
In |i~5[| , [17 1, Kiihn combines restricted annealing with the phenomenological renormaliza- 
tion group method to study the critical behaviour of the 2-D spin-diluted Ising model, 
while Kiihn et al. [|l(| apply constrained annealing to provide an exact solution of van 
Hemmen's spin-glass model [|l~8|]. 

Somewhat earlier restricted annealing was rediscovered by Schwartz [19] and put for- 
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ward in the debate about the lower critical dimension of the random field Ising model 
(RFIM), incidentally producing results compatible with the more recent rigurous proof of 
stability of ferromagnetism in the three dimensional RFIM [2G|. The most recent rediscov- 
ery of the method is due to Serva and Paladin |^] , who used it to study the mean-field 
version of the RFIM, which is actually somewhat simpler that the model considered in 
[16], and who produced approximate numerical solutions of the RFIM in 1-D. 

The purpose of the present paper is to give an outline of the formal and systematic 
aspects of Morita's equilibrium ensemble approach to disordered systems. Sec. 2 summa- 
rizes for future reference the general theory, and formulates constrained annealing as the 
canonical approximative scheme that derives from it. In this context it is emphasized that 
the equilibrium ensemble approach is nothing but a generalization of the grand ensem- 
ble idea jT^]. Alternatively, it can be looked upon as a method for obtaining variational 
lower bounds for quenched free energies. In Sec. 3 we specialize to bimodal bond- or site- 
disorder, and point out a connection to a different well established set of ideas, namely 
the relation to conventional low-concentration expansions and to perturbation expansions 
about homogeneous reference systems. For systems with bond- or site-dilution additional 
information about the exact disorder potential 4>(k) can be obtained which is used among 
other things to provide exact general solutions to the 1-D case. The equilibrium ensemble 
approach is well suited to study systems with correlated quenched disorder. Sec. 4 develops 
this point to some extent. In Sec. 5 we briefly indicate a range of alternative possibilities 
to formulate truncation and approximation schemes within the general approach, and we 
end with an outlook and concluding remarks. Constrained annealing as a method to study 
the thermodynamics and critical behaviour of disordered Ising models will be taken up in 
greater detail in an accompanying paper p2| , to be referred to as II in what follows. A 
short version collecting the main results of that study has appeared elsewhere [17]. The 
contents of Sec. 2 are not entirely new. Large parts of the general theory, in particular 
Sees. 2.1-2.4, have more or less explicitly appeared before — notably in ||, 113]. What we 
can perhaps claim originality for, is for pointing out heretofore unnoticed connections to 
other sets of ideas (to grand canonical ensembles, to variational bounding of free energies, 
and to low-concentration and perturbation expansions) as well as for the collection of 
exact results presented in Sees. 3 and 4. 



2 General Theory 

2.1 Equilibrium Ensemble Approach to Disordered Systems 

In order to avoid the complications associated with the evaluation of quenched averages, 
Morita Q suggested to introduce an equivalent ensemle in which configurational averaging 
becomes part of the thermal averaging procedure in an enlarged phase space that is defined 
over both, the configurational degrees of freedom k, and the dynamical variables proper. 
This is achieved by introducing a "ficticious" disorder potential cp(n) which is added to 
the system Hamiltonian 

H((t\k) ^H^(a,K) = H(ct\k) + <t>(K) (2) 
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and chosen such that the system described by Ti^ will have thermodynamic equilibrium 
properties identical to the nonequilibrium properties of the quenched system. Q In a first 
step, therefore, it is required that the Gibbs-Boltzmann distribution generated by TC^ 

p+{a,K) = ±exp[-pH*(<T,K)] (3) 
be equal to the non-equilibrium distribution 

P q (a,K) = ^re X p[-0H(a\K)] (4) 

describing the joint statistics of configurational and dynamical degrees of freedom in an 
ensemble of quenched systems. That is, one demands 

p^(a,K) =p q (a,K) (5) 

for all (a, k) . This equation can be solved for (j> to yield 

f3(f)( K ) = - ln[q(K)/Z( K )} -\d.Z* . (6) 

We shall in what follows refer to (|6|) as to the Morita equation. It expresses the fact that an 
equivalent equilibrium ensemble exists. As long as one is interested only in its probabilistic 
content, Eqs. (g), @ and (||) are all one needs in principle. In practice, of course, not much 
has been gained so far. The translation of the problem of configuration averaging in systems 
with quenched randomness into the language of an equivalent equilibrium ensemble is 
purely formal. It does not so far help to circumvent any of the difficulties for which it was 
invented to begin with. To use the idea in practice, one has to find a representation for 
the potential 4>(k), which is adapted to the given disorder problem, and which can be the 
starting point for the construction of systematic approximation schemes. 

Before turning to this topic in the following subsection, let us address the thermody- 



namic content of the equivalent equilibrium ensemble. It was pointed out by Huber [14], 
that due to the invariance of @ under the transformation 

0(«) -» 0(«) + <p , (7) 
Z^> Z*exp[-/J0o] , (8) 

the disorder potential 4> and thus the equilibrium ensemle's partition function Z^ are spec- 
ified only up to arbitrary constants <p$ (constant in the sense that they do not depend on 
the disorder configuration k) , which may, however, still depend on temperature, external 
fields, coupling constants, impurity concentration and so on. For the probabilistic content 
of the theory, this need not be a point of concern. If, however, one also demands that 

F* = -pr l \ogZ* (9) 

be identified with the free energy of the equivalent equlibrium ensemble, a normalization 
is clearly needed. It was demonstrated by Huber KM that thermodynamic consistency - 



1 Here and in what follows, we suppress indices signifying the system-size dependences of various quan- 
tities, in order to simplify notation. 
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equality of probabilistic and thermodynamic definitions of internal energy, entropy, and 
other thermodynamic functions involving derivatives of up to first order in temperature 
and fields — requires the configuration average of <f> to vanish, 

(^)) ? = E?w^) = ' ( 10 ) 

K 

Indeed, with this condition it follows from @ that 

-^ = lnZ* = (InZ(«)> g -(ln g («)>, , (11) 

so that gives the quenched free energy plus a contribution which is easily identified 
with an entropy of mixing, 

kg 1 So = - (lng(«)>, = -5>(«)hg(K) . (12) 

K 

Since this extra contribution is independent of temperature and external fields, it is ir- 
relevant for thermodynamics. With the normalization condition imposed, the disorder 
potential <p reads 



(13) 



Z(k) \ Z(K), q 

which shows that a full specification of <ft requires the computation of the quenched free 
energy — precicely the quantity the computation of which the equilibrium ensemble ap- 
proach was invented to facilitate, not to speak of the difficulties of solving or analysing 
models with a contribution to the potential energy as complicated as (p|). On the other 
hand Eq. ( |l3|) shows that, having computed a properly normalized 4>(k), it suffices to 
evaluate it for a simple homogeneous configuration k to obtain essentially the Brout free 
energy. 



2.2 Constrained Entropy Maximization 

In order to introduce an alternative and more explicit representation of the disorder po- 
tential 4>, which can serve as a starting point for systematic approximation schemes, it is 
advantageous to recall Mazo's |2j| information-theoretic justification of Brout's averaging 
prescripton. 

The distribution p q {a,K) which reproduces Brout's averaging prescription for the free 
energy can according to Mazo be derived from Jaynes' information-theoretic approach to 
statistical mechanics PH - That is, it is given as the unique distribution which maximizes 
the entropy functional 

S[p] = -k B J2p(a, k) lnp(cr, k) (14) 

under the constraints 

(H(*\K)) p =E., K p(°,K)n(a\K) =E (15) 
(S(k-k')) p =£*,„?(*,«)<$(«-«') =g(«'), Vac'. (16) 
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Eq. (H) fixes the average energy of the system, while (jl6|) demands that the search for 
maximizing distributions is restricted among those compatible with the a-priori distri- 
bution q(n) characterizing the disorder. Introducing Lagrangian multipliers —k^fi and 
A;b(A(k) + 1) to enforce the constraints, one obtains p q (cr, k) defined by 

S\p q ] = max{S[p\; (n(a\ K )) p = E and (6(k - n')) p = ?(«') V «'} (17) 

in the form 

p q (a, k) = exp[-pH(<r, k) + A(k)] . (18) 
To satisfy dl6|), one on needs 

exp[A(„)] = f^- , (19) 
Z{k) 

where Z(k) is the partition sum at fixed disorder configuration k, so that ( pi ) does, indeed, 
agree with (f|), and we have seen in the previous subsection that it can be interpreted as 
an equilibrium distribution p^ in an enlarged phase space. 



2.3 Alternative Formulation of Constrained Entropy Maximization (Bi- 
modal Disorder) 

The alternative representation of <j> announced above is obtained by replacing the con- 
straints (fitf ) in the above entropy-maximization procedure by an equivalent set of con- 
straints, namely by the requirement that the search for maximizing distributions is re- 
stricted to those which reproduce the complete set of moments of q{n). 

To be specific, and to keep notation and formulae as trasparent as possible, let us first 
consider the case of systems with bimodal site (bond) disorder, in which a random variable 
associated with each site (bond) of the lattice can only take two values. We shall explain 
the necessary modifications to treat more general types of randomnes in Sec. 2.5 below. 
Bimodal disorder is conveniently described in terms of site (bond) occupation numbers 
k a = k a (n) which take the values or 1, if on the a-th site (bond) one or the other of 
the two possible realizations of the disorder is attained in the configuration k. The set 
of moments of q(n) is given by the expectation under g(/c) of all possible products of 
occupation numbers k a , 

U = (l[k a (K)\ , u,CA, (20) 

where u ranges through all subsets of the set A of sites (bonds) of the lattice, including 
the empty set. It is useful introduce occupation numbers for the sets u through 

k{k) = n k a {K) , (21) 

with k$(K) = 1, and to identify n with the set of vertices (bonds) for which k a = 1, that 
is, k = {k a ; a E A} = {a G A; k a = 1}. This allows to write 

f u = (*«(«)), = E ?(«) • ( 22 ) 
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In the present case, the equivalence of the system of moments f u and the distribution g(/c) 
follows directly from the invertibility of (p2|), 

(23) 



(*) = E(- 1 ) H " w /«. 



where \u\ and |/s| denote the size of the sets u> and k, respectively. The set of moments 
can be used to reformulate the constrained entropy maximization, by replacing (|l6|) with 



wCA. 



(24) 



Since (jig ) and (£4|) are equivalent sets of constraints, one can alternatively obtain p q (a,K) 
through 



S]p q ] = max {S\p\; {H{<j\k)) p = E and <A^(k)> p = f u , u C a} 



(25) 



Denoting the Lagrangian multipliers that enforce (15) and p4] ) by — /cb/3 and fcB(A tJ +5 lt , 0), 
one obtains the maximizing distribution p q in the form 



p q (a, k) = exp - /3H(ct\k) + ^ Xuk u (K) 



Writing 
one gets 



A© 



■InZ, 



(26) 

(27) 
(28) 



p q (a, k) = — exp -/3ft(er|«) + ^ A w ^(k) 

where the A^, a; / 0, must be determined such that p q satisfies (p4| ) for the set of moments 
of the distribution q(K) describing the disorder. Thus, an explicit representation of 4>(k) 
in terms of weighted sums of products of occupation numbers has been obtained. The 
representation (28) of the distribution p q (a,K) has been advanced already by Morita ||, 
in a more transparent way, though, by Sobotta and Wagner fTT 



Note that Z has the formal appearance of a generalized grand partition sum, with 
generalized chemical potentials given by fj, u = /? -1 A w . The corresponding generalized grand 
potential is 

n = -P~ 1 lnZ, (29) 



and one has 



(30) 



Returning to the general considerations of the previous subsection, we realize that in the 
form in which the disorder potential appears in (p8|), 



(31) 
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it does not satisfy the normalization condition (|10|). To enforce it, one changes <j) according 
to 

0<j>{n) -> - MM*) ~ U) • (32) 

On the level of thermodynamic potentials, this change amounts to a Legendre-transform 
from the generalized grand potential to the corresponding free energy, 

f)C) 

n^ F t = n + J2 /WL = n - E ^jr~ ■ ( 33 ) 



Finally, using (p7[), ( |33D and the fact that = 1 due to the normalization of = p g , 
one obtains an intriguingly simple representation of the equivalent ensemble's free energy, 
viz. 

-)3F^ = \nZ^ = -J2 <Wu, ■ (34) 



Note that determining the A w such as to make p^ satisfy the constraints defining the 
disorder distribution is tantamount to minimizing InZ^ with respect to the A w . Phrased 
differently, it means minimizing In Z<t> over a space of functions 4>(n) which are of the form 



2.4 A Canonical Approximative Scheme: Constrained Annealing (Bi- 
modal Disorder) 



The simplicity of ( p4[) is, of course, purely formal and hence deceptive. Nevertheless, some- 
thing non-trivial has been gained, namely an explicit representation for the disorder po- 
tential <j), together with a method to compute it that immediately lends itself to the 
formulation of a systematic scheme of approximations. This scheme consists of replacing 



the full set (24) of constraints by a subset which requires the maximizing distributions to 
reproduce only a given subset of the full set of moments of q(n). 

Formally, let Q be a subset of the powerset V(A) of the vertices or bonds of the lattice, 
and define p^ as the unique distribution that maximizes the entropy functional (|l4|) under 
the constraints 

(H(a\K)) p =Z«, k p(o;k)H(o-\k) =E (35) 

(M K )) P =Y,a, K p( a , K ) k U K ) = fuj , ^eQ, (36) 

that is, 

S[p Q ] = max{%]; {H(a\n)) p = E and <M«)> P = fu> , u E Q\ . (37) 

Introducing Lagrangian multipliers —k-^fi and &b(A5 + 5^$), u E Q, one obtains the 
maximizing distribution p® in the form 



p Q {a, k) = exp 



PH(v\k) + X; A£M*)] • (38) 



S 



Writing 

one observes that this takes the form of a generalized grand-canonical distribution 



\f = -lnZ Q , (39) 



p Q (a, k) = exp [ - PH(<t\k) + £ A£*u,(«)] , (40) 



where one has to determine the A® = such as to satisfy the constraints (^). In order 
to obtain the free energy corresponding to the Q-approximation, one has to enforce the 
zero-mean-potential condition ( |i~0| ) by changing <p^ according to 

/^(k)-- ^S(kM-U) ■ (41) 

As above, this amounts to a Legendre transformation from the generalized grand potential 

n Q = -p~ 1 lnZ Q (42) 

to the corresponding free energy 

Once more, we have a simple representation for the free energy, 

-PF Q = \nZ Q = -Y, A£/ w . (44) 



Again, to determine the A® so as to make satisfy the restricted set (|36|) of constraints 
is equivalent to minimizing In with respect to these parameters or, put differently, over 
a space of functions 4>(n) which are of the form (|4l"|). 

Unlike before, the simplicity of (p4|) need not be purely formal. For instance, if one 
restricts one's attention to simple approximations which reproduces only a small subset 
of the set of moments of q{n) — say the translationally invariant expectations (k a ) and 
the translationally invariant nearest neighbour correlations (k a kb) in the case of bimodal 
site-disorder, so Q = {{a} ag A> {(a, fr)}a,&eA}> with (a, b) denoting nearest neighbours - 
then the corresponding Hamiltonian T$ acquires only on-site terms and nearest neighbour 
interactions between the k. So one obtains a fairly simple translationally invariant model. 
For its study the full machinery of equilibrium statistical mechanics and an established 
arsenal of approximation methods is, at least in principle, available. We shall demonstrate 
below that such simple approximations can already yield very precise, if not exact, results. 

There are a number of requirements the set Q should meet to give rise to a sensible 
approximation for the underlying problem of describing systems with quenched disorder. 
First, to ensure that the distribution is normalized, Q must contain the empty set, 
since the normalization constraint is expressed precisely by the condition /g = {k$) = 1. 
Moreover, if the disorder is homogeneous in the sense that the distribution q{n) has a set 
of moments invariant under the group Q of operations which map the underlying lattice, 
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hence the set A of sites (bonds), onto itself, the set Q should be such that it respects these 
symmetries, i.e., it ought to be closed under the operations g of the symmetry group Q 
of the underlying lattice. Formally, if uj G Q then gui = {a G A; g~ 1 a G uj} should also 
be a member of Q. The invariance of the set of moments under Q entails a corresponding 
invariance of the Lagrangian multipliers, A^ = A^ for all g in Q. In other words, the set 
of constraints can be decomposed into equivalence classes under Q. 

Fixing Q so as to concentrate on moments of the disorder distribution which are 
of finite order or which involve occupation numbers referring to sites (bonds) only up to 
some maximal distance obviously entails a notion of locality for the corresponding disorder 
potential 0®(k). In what follows, we shall therefore refer to approximations of this type 
as Q-local. 

2.5 Exact Lower Bounds for the Quenched Free Energy 

Before turning to a more detailed evaluation of these ideas, let us not fail to mention that 
the approximative scheme described above provides exact lower bounds for the quenched 
free energy. 

This is easily seen by realizing that the constrained entropy maximization scheme 
described above can be substituted by a constrained minimization scheme for the free 
energy functional 

F[p] = Y^P^^n^^+p-^Pi^K) In p(a,K) (45) 

using the set of moments (k^) = / w , uj G Q as constraints. This yields 

F[p Q ] = mm {F\p\; <^(«)) p = f u uj G q| . (46) 

With suitable Lagrangian multipliers, the p Q so defined coincides with (pQ) as obtained 
from the entropy maximization scheme. If, moreover, one enforces the zero-mean-potential 
condition @, then F[p Q ] = F Q , with F Q given by (||). 

The dependece of on the set Q C T-'(A) is monotone in the sense that 

F Q < F Q! < p<p if q c Qf c V{K) (47) 

so that F® approaches the quenced free energy F^ monotonically from below, as the set 
of constraints, as characterized by Q is increased. A different way to see this is to recall 
that the lnZ s are obtained by minimizing In over a space of Qdocal functions 4> — a 
function space which becomes larger as the set Q is enlarged. 

2.6 Non— Bimodal and Continuous Disorder Distributions and the Vari- 
ational Content of the Equilibrium Ensemble Approach 

The modifications necessary to obtain explicit representations of the disorder potential 
and approximative schemes for non-bimodal or continuous disorder distributions q{n) can 
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be inferred from what has been outlined in Sees. 2.3-2.4 above. We will briefly state them 
in what follows. They also point — more obviously than our previous considerations - 
to the variational content of the equilibrium ensemble approach to disordered systems. 

For the remainder of this section, we drop the restriction that the k a be bimodal. 
To begin with, let us assume for simplicity that all moments of the (non-bimodal or 
continuous) disorder distribution q(n) exist. They are given as expectations under q(n) of 
all possible products of powers of the k a , i.e., by 



/"=( TT *£"(«)) ,^GA, n a GW, (48) 




with n = (n a ) aeuJ . The partition function of the equivalent equilibrium ensemble describing 
the disordered system is then given by 

= J2 exp[-/JW(<r|/«) - (3(/)(k)} , (49) 



with 6 of the form 



m*)= E ww-/:). (so) 



Here we have introduced 

TO = 11 *£"(«)■ ( 51 ) 



In (|49|), ( p0[) we have chosen to satisfy the zero-mean-potential condition (10) at the 
outset. The copling constants A" have to be determined such that the Gibbs distribution 
p4> corresponding to (|49| ) reproduces the moments of the underlying disorder distribution 
q, i.e., (k™) = f™ is supposed to hold for all to and n. This amounts to minimize lnZ^ 
as a function of the A", in complete analogy with the bimodal case. As above, the A" can 
be interpreted as Lagrangian multipliers of a constrained entropy maximiation scheme or, 
alternatively, of a constrained free energy minimization scheme. 

As before, approximations are obtained by enforcing only a subset of the full set of 
constraints, that is, by requiring the Gibbs distribution to reproduce only a subset of 
the full set of moments ( |4"8| ) of the disorder distribution q(n). This can be done in several 
ways, for instance by requiring a match with the f™ only for u in some subset Q of the 
powerset V(A) of the set of sites (or bonds) of the lattice, by matching moments only 
up to some maximal value for the n a , or by combining these two truncation schemes in 
various ways. 

If one defines the approximations through a moment matching scheme restricted only 



with respect to the choice of the oj in (|4q) and (50), demanding that oj G Q C V(A) as 
in Sec. 2.4, the minimization of lnZ^ with respect to the correspondingly restricted set 
of parameters A", lo G Q, amounts to a minimization over a true infinite dimensional 
function space, i.e., as a variational problem. It amounts to minimizing YnZ^, interpreted 
as a functional of <fi over a space of functions (j> which are Q-local in the sense defined 
above. That is, they are of the form 

<k«)= E (^(«)-(va«)> 9 ) , (52) 
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each ip u being a function of the k a , a G uj, which is analytic in every component, since 
power series expansions were assumed to exist; cf. Eq. (|50|), Restricting also the range 
of n a values would further restrict the set of functions to multinomials of some maximal 
degree. 

Conversely, dropping the analyticity constraint on the ip^ will enlarge the function 
space sufficiently to treat cases where not all moments of q(n) exist. 

Having said this much, it should have become clear that the equivalent equilibrium 
ensemble approach could have been formulated as a variational problem right away. To 
wit, consider Z^, interpreted as a functional of 

= Z[<f>] = £ Z{k) exp[-^(«)], (53) 

K 

and let us assume at the outset that the </> under consideration satisfy the zero-mean- 
potential condition (|Toj) ; alternatively, one might impose it a-posteriori. Using the disorder 
distribution q[n) this may be rewritten as an average: 

Z* = Z[<j>] = (Z(k) exp[-P<l>(K) - In </(«)]>, • (54) 
By Jensen's inequality, 

IrZ* = lnZ[cf>} > (In Z( K )) q - (lnq( K )}) g , (55) 

where we have invested {4>(n)) q = 0. The task then is to minimize the left hand side of 
(pa) over a suitable space of functions (ft satisfying (<P(k)) = 0. This variational point 



of view was utilized by George et al. ]10[ in their study devoted to the computation of 
phase-boundaries of bond-disordered 2-D Ising models — without connecting it to the 
equilibrium ensemble idea and to the heuristics associated with it, and without, so it 
seems, the awareness that it allows to proceed to exact solutions. 

Indeed, the inequality in (|55| ) can be satisfied as an equality, i.e., the lowest bound 
attained and thereby the exact quenched free energy computed, if (/)(k) is chosen such that 

Z{k) exp[— /3(P(k) — lnq(K)] = const (56) 

in the sense that this quantity does not depend on k. The value of this constant must, of 
course, be Z 1 * ', so (|56|) is nothing but the Morita equation @ in disguise. That is, it is the 
formal solution of the variational problem 

^-0. (57) 



Approximations are obtained by restricting the space of functions <j) over which a 
minimum of lnZ[</>] is sought, in other words, the domain of defenition of the functional 
In Z [</>], to functions which are Q-local in the sense explained before. Approximations can 
systematically be improved by increasing the size of the set Q 
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3 Recursive evaluation of coupling constants 



For systems with quenched, uncorrelated bimodal site or bond disorder, one can write down 
more or less explicit expressions for the coupling constants A^ of the full disorder potential 
4>. It turns out that these expressions generate conventional low-concentration expansions 
for the quenched free energy @) as well as perturbation expansions about pure reference 
systems. In the present section, we shall establish this connection. In Sec. 3.2, we specialize 
to models with bond- or site-dilution for which additional exact information about the 
exact disorder potential 4>{k) is available. Among other things, we use it to reproduce 
the known exact solutions of one-dimensional models with quenched site or bond dilution 
and to exhibit simplifying features of these solutions in the absence of symmetry breaking 
fields. 



3.1 General results 



We consider systems with quenched, uncorrelated bimodal site or bond disorder, defined 
on a lattice C containing N sites. As before, let A denote the set of sites (bonds) of the 
lattice. The distribution q{n) describing the disorder configurations 



K = {k a ; a G A} = {a G A; k a = 1} 
is characterized by a single parameter p, 

?(«) = J] p fe -(l - p) 1 ^ = pN(i _ p) |A|-|«| 
aeA 



so that 

We use the Morita equation 



fui — {ku)) q — P 1 '• 

AD, 

j3(f)( K ) = - hx[q(K)/Z{K)} -XnZ^ 



(58) 
(59) 
(60) 



(61) 
(62) 



and the representation 

- - fw) 

for the exaxt potential 4>(k), to obtain the coupling constants A^. This is achieved by 
inserting (^2|) and (|60|) into the Morita equation, and by considering this equation first for 
k = 0, which gives 



/ty(0)= Y, = -|A|ln(l-p)+lnZ(0)-ln^ , 



(63) 



and then for k ^ 0, which yields 

04>{k) 



am- E A - 

-|A|ln(l - p) - \n\hx 



1-p 



+ In Z(k) - In 



(64) 
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Subtracting (|6J) from (|63|) , we get 



|k| In 



1 



lnZ(«) +lnZ( 



(65) 



To determine the A w individually, we introduce the quantities x u , uo C A, through 

^c«) = n x - ( 6e ) 



or, equivalently 



InZ(Ae) = ^ In 



(67) 



where k ranges through all subsets of A (including the empty set) and the product in (66) 
is over all subsets of k, including again the empty set. Eq. (^) may be interpreted as a 
recursive definition of the alternatively, one may invert this relation directly to yield 



lnx w = ^(-)l w H K llnZ(K) 



(68) 



Inserting this into (p5|), and recalling that is supposed to hold for all k ^ 0, we can 
immediately conclude 



A, 



In 



1 



'M,i 



(69) 



This is the more or less explicit expression for the coupling constants A w of the exact 
disorder potential <j) announced above. Using this result, and returning, once more, to 
the k = version ((6^) of the Morita equation, we get an (again more or less explicit) 
expression for the free energy of the equivalent equilibrium ensemble, 



hiZ+ = -|A|[plnp+ (1 - p) ln(l -p)]+J2 <° M lnx - 



(70) 



uCA 



The first contribution in (|70| ) is readily identified as the expected contribution of the 
entropy of mixing. Using (pq), the second contribution can be cast into the form 



£ pH In xu= ^pl K l(l-p)l A H«linZ(«) . 



(71) 



jCA 



kCA 



Up to the factor — this is just the Brout avarage (|l]) of the free energy for the case of 
quenched, uncorrelated bimodal site or bond disorder as it should, since no approximations 
were involved in the above considerations. 



Clearly, we cannot expect to be able to evaluate (|70|) completely, if we are unable to do 
so in the more standard formulation (|ll]). Nevertheless, we have gained something useful. 

First, Eq. (|70|) together with (^) amounts to a formal model-independent prescription 
for reorganizing the Brout average in such a way that it is turned into an expansion ordered 
by increasing powers of p. That is, for systems with bimodal bond- or site-disorder, 
there is apparently an intimate relation between the equilibrium ensemble approach to 
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disordered systems and conventional low concentration expansions for such systems. It 
should be noted that we were as yet not forced to specify the physical meaning of (low) 
concentration in the above considerations. Clearly, this will depend on which meaning we 
decide to associate with the k a = 1 states of the disordered system. Think for example of 
a randomly spin-diluted magnet. If we choose k a = 1 to designate that site a is occupied 
by a spin, then p has the meaning of a (spin) density, and (|7C|) gives indeed a conventional 
low concentration expansion. We might, of course, also have chosen k a = 1 to designate an 
empty site (occupation by a vacancy), thereby generating an expansion about the pure, 
homogeneous reference system. Clearly these different choices are giving entirely different 
meanings to the expressions appearing in (p8|) , (70) and (|7l|). The quantity = lnZ(0), 
for instance, gives the free energy of the empty lattice, if k a = 1 denotes occupancy by a 
spin, whereas it gives the free energy of the fully occupied homogeneous system, if k a = 1 
denotes that site a is empty. 

Second, we can use the formal results obtained above to compute at least the simplest 
coupling constants of the exact disorder potential explicitly. This turns out to provide 
useful information for the investigation of the approximative systems introduced in the 
previous subsection. 

Third, for models with bond- or site-dilution, we will be able to show in Sec. 3.2 
below, that a large nontrivial class of coupling constants of the exact disorder potential <f> 
will vanish in the absence of symmetry-breaking fields. This piece of information, too, can 
be useful as a guide for finding efficient sequences of approximations within the scheme 
presented in Sec. 2.4. 



3.2 Models with Bond- and Site-Dilution 



Models with bond- and site-dilution exhibit a simplifying feature which allows to carry 
the formal developement of the equilibrium ensemble approach one step further, namely 
the notion of independent non-interacting clusters in terms of which every disorder con- 
figuration can be characterized. 

If we adopt the convention that k a = 1 denotes an occupied site (bond), so that the 
equilibrium ensemble approach generates a conventional low concentration expansion as 
explained in the previous subsection, we can establish the following results: 



(i) x^ = 1, hence A w = 0, if co denotes a disconnected set of occupied sites (bonds). 

(ii) In the absence of symmetry breaking fields, x wU { a } = 1, hence \ w u{ a y = 0, if uj denotes 

a connected cluster containing at least two occupied sites (one occupied bond) to 
which a single site (bond) a is added to form a dangling connection with to. 

(iii) In the absence of symmetry breaking fields, x^ Y . W2 = 1, hence A^.^ = 0, if uj\ ■ UJ2 
denotes a connected cluster of sites (bonds) which consists of two separate sub- 
clusters oj\ and L02 that are conneted only via a single common site. 



Note that (ii) can in fact be read as a special case of (iii). All three results hold for 
site or bond diluted models with nearest neighbour interactions on arbitrary graphs. They 
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follow directly from the definition (]66| ) of the x w . We establish them separately for systems 
with site-dilution and with bond-dilution. 

3.2.1 Models with site dilution 

Identifying the configurations k of the disordered system with the sets of occupied sites, 
we have 

Z(0) =x 9 = l, (72) 

and 

Z(k)= J] x u (73) 

for k ^ 0. The proof of (i) is by induction. First, using Z({a}) = Xf a x = Z\ for configu- 
rations in which only a single site a is occupied, and Z({a} U {b}) = Z({a})Z({b}) = Z 2 , 
if a and b are not nearest neighbours, one immedieately obtains £{ a ,6} = 1 f° r the sim- 
plest disconnected set {a, b} containing only two sites. Supposing now that (i) has been 
proven for all disconnected sets u containing up to n sites, and assuming that k denotes 
a disconnected site containig n + 1 sites, one has 

Z(k) = J] Z( Ka ) , (74) 

where the product is over all independent clusters (maximally connected subsets) of k. By 
supposition 

z(k) = n z m = n ( n 

LOa) x K (75) 

entailing x K = 1 also for the disconnected set k containing n + 1 sites, which proves (i) for 
the case of site-dilution. 

Next, to prove (ii) we consider configurations of the form k = wU {a}, where the site a 
forms a dangling connection with uj. Denote by E(a a ,ab) the interaction energy between 
the degree of freedom sitting on the dangling site a and erj,, the degree of freedom in uo to 
which it is connected, in the absence of symmetry breaking fields. This interaction energy 
is assumed to be invariant under an internal symmetry group of the model in the sense 
that 

i = 5>xp[-/3£(a a ,cT 6 )] (76) 

is independent of at, so that the partition function Z(k) = Z(u U {a}) in the absence of 
symmetry breaking fields is given by 

Z(k) = Z{u U {a}) = Z Z(u) . (77) 

The proof of (ii) is likewise by induction. First, we have Z({a}) = Z\ = X{ a y and, by 
the factorizing property fl77D, Z({a,b}) = Z2 = Z Z\ = S{ a ,6} f° r a configuration 
containing only two neighbouring sites a and b. This gives ^{ a ,6} = Z jZ\. Adding another 
dangling connection in k = {a,b,c} gives Z({a, b, c}) = Z3 = Z 2 Z\ = x^ a y x 2 ab y <E{ ai & )C } ; 
so X{ a 5 c \ = 1, which verifies (ii) for the smallest conceivable set of the form oj U {a} with 
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ijj containing two sites. Suppose now that (ii) has been established for all sets of the form 
ijj U {a} with a dangling connection from tu to a single site a and with uj containing up to 
n sites, and let k be of the form k = u U {a} with oj containing n + 1 sites. Then 

Z(k) = Z(u U {a}) = Z Z{u) 

= ( II ( II X ^'U{a}) ^{a} Xu,U{a} > (78) 

SO 

(79) 

By supposition all avuja} in the above product except £{ a ,6} are unity, where b denotes 
the unique site in Co connected to a via the dangling connection. Since X{ a b y X{ a } = Z, we 
can conclude that ££u{a} = 1 as well, which completes the proof of (ii) for the site-diluted 
system. 

To verify (iii), note that the result X{ a ^, c } = 1 f° r three neighbouring sites o, b, c 
obtained above establishes (iii) for the smallest conceivable set of the form lo\ ■ U2, with 
uj\ = {a,b} and u>2 = {b, c}. The result follows by induction after noting that in the 
absence of symmetry breaking fields Z{ui\ -002) = Z\ Z (oji) Z (002) , and by considering the 
consequences of this identity in the ^-representation. 



3.2.2 Models with bond dilution 

In the bond-diluted case, we identify k with the set of occupied bonds. Thus, if k = 0, 
the system consists of a collection of N non-interacting degrees of freedom, and hence 

Z(0) = Zf , (80) 

where, as before, Z\ denotes the parition function of an isolated degree of freedom. If k ^ 0, 
there will be a set of interacting degrees of freedom — those sitting in clusters formed 
by bonds occupied in At — but in general also a collection of non-interacting degrees 
of freedom. Denoting by Z K the partition function of the degrees of freedom interacting 
through the set of bonds occupied in k, we can write 

Z{k) = Z» ss < ( 81 ) 

where n(n) is the number of vertices linked by bonds occupied in k. Returning to the x u 
representation ([66|), we identify 

£0 = (82) 

and 

Z^ = ^k= II *«■ (83) 

The Z(k) have properties analogous to the Z(k) in the site-diluted problem. In particular, 
if k = K\ UK2 and k\ and K2 are mutually disconnected sets of occupied bonds, then Z(k) = 
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Z{k\)Z{i^i). This factorizing property again entails x^ = 1, if u denotes a disconnected 
set of bonds. The proof is by induction in complete analogy to the site diluted case, and 
it will not be repeated here. 

Next, to establish (ii) for the bond diluted case, we note that 

Z(lo U {a}) = — Hrrr = ~ rr = tt^US) (84) 

in the absence of symmetry breaking fields, if a denotes a dangling connection added to 
the set lo; here Z has the same meaning as in the site diluted case. Again, by inductive 
reasoning in complete analogy with the site diluted problem, one verifies (ii), i.e., i wU | a } = 
1 if a; contains at least one occupied bond. 

Finally, to prove (hi), we note that 

Z(wi • wa) = , ^7 , = = Zi °,\ Zu ;\ = Z{u x )Z{u 2 ) , (85) 

where the second equality requires the absence of symmetry breaking fields. This factor- 
izing property of the Z{uj\ ■ 102) entails (hi) on the level of the x w -representation. 

3.2.3 One dimensional models with bond or site dilution 

One dimensional systems with bond or site dilution are particularly simple because in 
1-d one has a complete overview over the family of connected clusters which consists just 
of isolated (finite) chains. This allows to give an explicit representaion of the x u defined 
through (|6^) in terms of chain free energies, which we briefly record here. 

For the site diluted problem, denoting by Z n the partition function of an isolated open 
chain of n occupied sites, and by x n the corresponding x variable, we have 

Zn = ft x n k +1 ~ k , (86) 
fc=l 

which gives x\ = Z\, and (with Zq = 1) 

x n = Z n Z n _ 2 /Zl-i , n>2. (87) 



Inserting this into the expression ( |70|) for the free energy of the equivalent equilibrium 
ensemble and normalizing with respect to sample size N = |A|, we get 

-Pf* = lim jV -1 ln^ 

N^oo 

00 

= -plnp-(l-p)ln(l-p) + (l-p) 2 Y / P n ^Z n . (88) 

n=l 

Up to the contribution coming from the entropy of mixing, this reproduces the exact 
solution as obtained, e.g., by Wortis [25]. In the absence of symmetry breaking fields we 



have lnx n = for n > 3, so in this case 

= -plnp-(l- / o)ln(l-p)+plnZi+ / 9 2 ln(Z 2 /Z 1 2 ) . (89) 
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In the bond diluted problem, isolateds chain containing n successive bonds constitute 
chains of n + 1 interacting degrees of freedom. If we formulate (|83| ) for such chains, it 
reads 

7 n 

fs- = n > oo) 

z i fe=i 

where the Xk now refer to chains containing k bonds. This gives x\ = Z2/Z 2 , and 

X n = Zn+xZn^x/Zl , n > 1. (91) 



Comparing with the corresponding expression ( |87| ) for the site-diluted case, we note that 
x bond _ x site^ £ Qr n > 1 p r the free energy per site in the bond-diluted system we thus 
get (recalling x$ = Z^) 

00 

- /3f* = -plnp-(l-p) ln(l - p) + (1 - p) 2 P n ~ X ^Z n , (92) 

71=1 

thereby reproducing again the exact solution^]. In the absence of symmmetry breaking 
fields, lnx n = for v > 2, and the expression simplifies to 

-/3/^ = - / 9lnp-(l- / 9)ln(l- /0 )+lnZi + pln(Z 2 /Z 1 2 ) . (93) 

Eqs. ( p8|) and (|92] ) can be evaluated for all models for which the finite chain thermody- 
namics is known, e.g., from transfer-matrix techniques. 

Note that the result (ii) concerning the vanishing of coupling constants of the disorder 
potential of the form X w ij^ a y in the absence of symmetry breaking fields implies that in 
this limit very simple approximations in the moment-matching approximation scheme 
described in Sec. 2.2 already provide exact solutions for bond- or site-diluted chains. For 
the bond-diluted case, it is the simplest conceivable annealed approximation enforcing 
only the average occupancy (k a ) = p of the bonds. In the case of site-dilution, one has to 
fix the average site-occupancy at p and the expectations (k a kb) for neighbouring sites at 
p 2 in order to reproduce the thermodynamics of the quenched system exactly. 

On the contrary, if there is a nonvanishing symmetry breaking field, then Eqs. (87) 



and ( |9l[) imply that no finite approximation in the moment-matching scheme described 
in Sec. 2.2 is exact. Nevertheless, the simplest approximations enforcing only (k a ) = p in 
the bond-diluted case, and both, the site occupancy {k a } = p and the nearest neighbour 
expectation (k a k a+ i) = p 2 in the site-diluted case already give rather precise results, as 
will be demonstrated for the Ising model in II. 



4 Thermally Correlated Frozen-In Disorder 

Having treated uncorrelatded bimodal bond- and site-disorder in some detail, it is perhaps 
appropriate to recall that the equlibrium ensemble approach to disordered systems is not 
restricted to this case. Except in the previous section we have, in fact, never made use 
of the specific simplifications pertaining to uncorrelated disorder. The moment matching 
schemes explaind in Sees. 2.4-2.6. can naturally also be put to work in the case of correlated 
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disorder. The only nontrivial piece of information needed is contained in the system of 



moments ( |20|) or (48) for the bimodal and non-bimodal cases, respectively. In real physical 
systems, which do as a rule exhibit some degree of correlations between the configurational 
degrees of freedom, a subset of the full set of moments, obtained e.g. through scattering 
experiments, is often the only piece of information that is actually available about the 
disordered system in question. So, in situations where relatively simple approximations in 
a moment matching scheme can be expected to produce reliable results, the equilibrium 
ensemble approach appears to be an ideally suited tool to use. 

In the present section we will treat correlated bimodal bond- and site-disorder in some 
detail. We will restrict our attention to the case where the correlations are produced as 
a result of some previous annealing proccess (in the narrow metallurgical sense). That 
is, we suppose that the quenched disorder distribution is an equilibrium distribution pro- 
duced at values of the external parameters (temperature, fields, chemical potential, etc.) 
characteristic of the preparation process, but in general different from the values at which 
experiments on the system are actually being perfomed. 

Let us assume that the system is prepared at (inverse) temperature f3 and that it is 
described by the total Hamiltonian 

H tot (o-,K)=H(o-\K) + V(K) (94) 

where 7i(a\K,) is the Hamiltonian of the disordered system describing the energy of the 
cr-degrees of freedom at fixed disorder configuration k, albeit at parameter settings char- 
acteristic of the preparation process, which we indicate by the overbar. The potential V(n) 
accounts for the fact that there may be an additional interaction energy between the con- 
figurational degrees of freedom which produces correlations between the k a over and above 
that mediated by a. Think, for example, of a dilute magnetic alloy. The distribution of 
particles will usully be governed, both, by chemical interactions — accounted for by V(k) 
- and by K-dependent interactions between magnetic degrees of freedom described by 
H(o"|k). In typical cases, the former may, in fact, be expected to dominate. 

We assume that during the preparation process the temperature is high enough to 
allow the system to come to complete equilibrium. The distribution of the configurational 
degrees of freedom is then given by 

?(«) = ^exp[- / 9y(«)] , (95) 

where Z{k) = J2a ex P[ — /37~C(o~\k)]. If the system is then cooled down to lower tempera- 
tures, the impurity distribution fl95|) may be regarded as fixed, and the disordered system 
characterized by this quenched distribution will exhibit frozen correlations between im- 
purities. In what follows, we will show that the formal developments presented in Sec. 3 
for the uncorrelated case can be carried almost as far in the case of thermally correlated 
quenched disorder. 

To this end we return to the Morita equation (|6|) and the representation 

/ty(«) = - ]T XUkM-f w ) (96) 
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of the disorder potential <ft(n), the overbar indicating that the moments are evaluated with 

The coupling constants of the disorder potential for this problem of correlated disorder 
can now be determined along the same lines as in Sec. 3.1. For the bimodal disorder type 
considered here, we can always write 

V(k) = £ vMk) (97) 

Using the x representation (|6G|) for both Z(k) and for the partition function Z(k) with 
the parameters settings prevalent during the preparation process, we get 

Ao, = —ftv u - hi(x u /xu) , oj^% . (98) 

Returning to the k = 0-version of the Morita equation, and recalling that = 1 irrespec- 
tively of the correlations, and that Z($) = x$ and similarly Z{%) = x$, the free energy 
is obtained in the form 

-0F* = In = In Z + ]T &>u> + MW^))/« 

uCA 

= -(lnq( K )) q + ^hixj* . (99) 

ujCA 

The reader is invited to compare this with @. The main additional obstacle on the 
way to evaluating the moment expansion ( p9| ) to high orders — assuming such an eval- 
uation were possible for the concentration expansion (|70|) in the uncorrelated case — lies 
in the which cannot be expected to be of such a simple form as in the case of un- 
correlated disorder. Nevertheless, analytic approximations are conceivable in the form of 
high-temperature expansions of the f u , if (3 is sufficiently small. In any case, as will be 
shown in II, relatively simple moment matching schemes, that is, simple approximations 
in the restricted annealing approach, which do not require knowledge of high moments, 
may still be feasible and useful. 

It should be noted that the disorder potential 4>{k) reduces to V(k), if parameter 
values are the same as those during system preparation, so that x^ = x u . The reader 
may convince herself that this just what is needed produce the correct limiting form of 
Z* in the case where q(n) is the equilibrium distribution of the configurational degrees of 
freedom k. 

Specializing to systems with bond- or site-dilution, we observe that the results con- 
cerning the vanishing of certain coupling constants of the disorder potential <f> stated and 
proved in Sec. 3.2 have simple analogues in the case of thermally correlated disorder. To 
see this, note that the results (i)-(iii), if formulated on the level of the x w , do not depend 
on disorder statistics or parameter settings; so they hold for the x w alike. For thermally 
correlated quenched disorder, statements to the end that x w = x w = 1 then translate into 
^uj = —P v uj- This implies a vanishing of the coupling constants as in the uncorrelated case, 
if the the potential V(k) involves no interactions beyond those coupling configurational 
degrees of freedom on neighbouring sites (bonds). 
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Taking site-diluted chains in the absence of symmetry breaking fields as an example, we 
obtain the following surprisingly simple expression for the free energy per site as computed 
in the equilibrium ensemble approach, 

- = k^s + f x In Zx + h HZ 2 /Zf) , (100) 

where sq is the entropy of mixing per site, while fx = {k a )q = p and fi = (k a k a +x)q- That 
is all effects of quenched correlations on the free energy of the system make themselves 
felt only through the modified nearest neighbour correlation. This ceases to be the case 
if a symmetry breaking field is applied. Note that symmetry breaking fields during the 
preparation process manifest themselves in this context — beyond their effect in modifying 
fx and /2 — only in modifying the entropy of mixing. 

Similarly, in the bond-diluted system with thermally correlated quenched bond occu- 
pancy one obtains 

- = fci^o + In Zx + fx HZ 2 /Zf) . (101) 

So the result depends only on the average occupancy fx = p of the bonds and is independent 
of correlations between them. 



The simplicity of the zero-field solutions ( 100 ) and ( |101| ) for the 1-D systems with 



quenched thermal correlations is not easily detected in the exact solution 26], which has 
been obtained by the direct averaging approach. 



5 Concluding Remarks 

We have presented an outline of the formal and systematic aspects of Morita's equilibrium 
ensemble approach to systems with quenched disorder ||. Hitherto unnoticed relations to 
other, more conventional sets of ideas, such as low concentration expansions, perturba- 
tion expansions about pure reference systems, variational bounding of free energies, and 
generalizations of the idea of grand ensembles have been pointed out and elaborated. 

The canonical approximation scheme within Morita's approach — a moment matching 
procedure that goes under the name of restricted or constrained annealing — is of a 
non-perturbative nature, which should also be clear in view of its relation to variational 
methods. Conversely, the expansion of the exact disorder potential was demonstrated 
to generate conventional perturpation expansions. A third, in some sense intermediate 
approach would be to attempt solving or analysing models with a truncated expansion of 
the full disorder potential exactly. 

For models with bond- or site-dilution, it was demonstrated that a large non-trivial 
class of coupling constants of the exact disorder potential </>(k) vanishes in the absence of 
symmetry-breaking fields. This fact may help to explain why relatively simple approxi- 
mations within the moment matching schemes or the variational formulation presented in 
Sec. 2 have been so successful in describing aspects of critical behaviour of e.g. disordered 
Ising models — undisputably so at least regarding the computation of phase transition 
lines (|| , ||, [10 1, [17]), still under debate though (see, e.g., p7j|), regarding the much 



more controversial and difficult question of critical exponents 117]] , 125 



22 



In the case of continuous disorder distributions, straightforward moment matching 
schemes as explained in Sec. 2 might not always prove to be the most efficient way of 
formulating approximative solutions. Instead of fixing moments up to a certain order 
directly, one might, for instance think of reproducing expectations only of certain com- 
binations of moments of the quenched disorder distribution in the equilibrium ensemble. 
The choice of particular combinations might be guided by attempts to exploit powerful an- 
alytic structures, such as provided by orthogonal function systems, or by physical insight. 
In bond-disordered Ising models, for instance, the representation of the exact solution in 
terms of weighted sums over configurations of van der Waerden polygons suggests that the 
decisive quantities to be reproduced in an equilibrium ensemble approach would be gauge- 
invariant loop-correlations of the natural high-temperature variables ty = tanh(/3Jjj), 
expressed here in terms of the conventional notation for the exchange couplings. That 
is, the disorder potential (f>(n) = 4>({tij}) would be designed to fix correlations of the t{j 

aroud closed loops C, i.e., correlations of the form (il(ij)cC a * their corresponding 
quenched values, each such correlation involving moments of the corresponding quenched 
{ Jij}-distribution of arbitrarily high order. This idea has been exploited by George et al. 
|jl(| to locate phase boundaries in disordered Ising models, and it is also being dicussed in 



a recent preprint of Paladin et al. [28|. 

Finally, it is perhaps worth pointing out once more that the equlibrium ensemble ap- 
proach appears to be well suited to study problems with correlated disorder, as is borne 
out by the general theory of Sec. 2, and by the more specific considerations of Sec. 4. 
Within the framework of (low order) moment matching schemes at least, correlated dis- 
order appears to be hardly more complicated than uncorrelated disorder. Regarding more 
conventional approaches, this can — at best — be claimed for correlated disorder of the 
Gaussian type. 
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